六类高超声速激波-激波干扰的数值模拟研究

基于粘性全N-S方程，利用有限体积法进行数值离散，对Ⅰ到Ⅵ六种类型的激波一激波干扰现象进行了数值模拟，其中采用了二阶Harten-TVD格式，LU隐式算法以及Baldwin-Lomax代数湍流模型。每种类型都得到了清晰的干扰流场结构，其中Ⅳ型干扰的模拟结果与实验结果进行了比较，符合较好。数值模拟结果的对比指出，由入射位置差异引起的不同类型激波-激波干扰结构有很大不同，对壁面压强分布的影响也存在较大差异。     

梯度加权局部分析解差分格式及其应用

提出了一种新的针对Ｎａｖｉｅｒ－Ｓｔｏｋｅｓ方程符合含有梯度项的对流扩散方程的梯度加权局部分析解的有限差分格式，且给出了通用的系数表达式。ＧＷ－ＬＡＤ格式是基于积分型的有限体积法并且通过附加源项的方法来实现。能够适用于粘性和无粘性及不可压缩和可压缩的压力驱动流动。通过对叶栅流道中二维湍流流的计算表明ＧＷＬＡＤ格式与ＥＤ格式相比，具有较高的精度和很快的收敛速度并且得到更好的压力场和速度场。     

涡喷发动机可靠性增长试验方案的设计

根据对产品可靠性的最低可接收值与幂律模型的可行性增长管理公式、可行性增长管理策略，推导了涡喷发动机可靠性增长试验所需的台数公式，参照某些涡喷发动机的试验数据及上述公式，给出了某B－2型涡喷发动机的可靠性增长试验方案及其理想增长曲线。     

高超声速流动分区对接网格算法研究

本文发展了一种能适应强间断流场的满足通量守恒的分区对接网格算法，并给出了保证通量守恒的通量插值系数应满足的条件。采用有限差分NND格式，求解层流NS方程，数值模拟了马赫数6—20绕圆柱的高超声速粘性流场。数值试验表明：方法具有良好的稳定性和收敛性，强激波能自由地穿越内边界。     

多重网格技术在复杂粘性流场计算中的应用及研究

流体力学控制方程的数值求解过程中,当网格加密、粘性效应强或来流马赫数低的时候,流场收敛非常缓慢。本文在多块网格、LUSGS隐式格式基础上采用多重网格技术,计算了NLR7301case#1两段翼型的低速绕流及DLRF4跨声速流场,在复杂二维、三维粘性流场求解中获得了满意的加速收敛效果。文中还比较了V循环、W循环的加速效率,为多重网格的工程应用奠下基础。     

三维自适应非结构网格的Euler方程解

 将 Ausm+ 迎风格式应用于三维非结构网格中求解 Euler方程。对单元变量进行重构以获得空间高阶精度,对时间域采用多步龙格库塔法推进,并采用了当地时间步长和隐式残差光顺技术来加速收敛。采用多点择优推进阵面法生成复杂曲面的三角形网格,利用推进阵面法生成四面体网格。采用网格自适应技术对网格进行局部加密,以减少总体网格数目,从而提高计算效率。最后给出了绕 ONERA M6机翼的跨音速流动及绕麻雀 导弹的超音速流动算例,结果表明了本方法的有效性。     

等价差分格式在三维复杂流场数值模拟中的应用

结合数值Riemann不变量和通量限制器（FluxLimiter）的概念，将基于Chakravarthy和Osher半离散差分格式的等价差分格式应用于三维复杂流场计算，应用时间相关法和全流场捕捉技术，获得了Ma＝10，a＝10°状态下球头绕流结果，壁面压力分布与Rusanov的计算结果相符合，在此基础上细致地刻画了在Ma＝7，a＝20°状态下绕流带座能飞船的复杂流场。     

驱动方腔内涡流的数值模拟

本文介绍一种求解涡度流函数形式的定常不可压缩Navier-Stokes方程的半隐式指数型差分格式。涡度方程用上述格式求解,而流函数方程用多层网格法求解。这种组合求解方式具有很好的稳定性及较快的收敛速度。本文对Re=100,400,1000,3000,5000的驱动方腔及Re=400,1000的驱动长方腔进行了数值模拟,结果与现有的计算相吻合。     

超声速主流中横向喷流场的激波—旋涡结构的数值模拟

本文利用NND格式,通过求解NS方程,对二维超声速主流中横向喷流干扰流场进行了数值模拟,计算清楚地给出了激波结构、回流区和混合层。本文计算得到的激波结构和实验相当一致。最有兴趣的是由于喷流的干扰,主流在喷口前发生主涡分叉,观察到三个流向旋转涡和两个反流向旋转涡;在喷口后的背风区,存在具有低压和回流区的尾迹。     

含运动边界二维多体干扰流场的数值模拟

本文采用二阶精度Godonov扩展方法求解非定常：Euler方程，以基于格心的自适应非结构网格的显式有限体积法为基础，空间和时间都是二阶精度，使用HLLC近似黎曼解的方法计算网格单元边界处的守恒量通量，对作俯仰振动的NACA0012翼型绕流及激波管中运动、静止圆柱的流场进行了数值模拟，并对数值模拟的结果进行了分析，为进一步对含多体分离流场进行数值模拟打下了基础。